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Abstract. The phase transition to mirrorless oscillation in resonantly enhanced 
O ' four-wave mixing in double- yl systems are studied analytically for the ideal case 

C_J I of infinite lifetimes of ground-state coherences. The stationary susceptibilities are 

obtained in all orders of the generated fields and analytic solutions of the coupled 
nonlinear differential equations for the field amplitudes are derived and discussed. 



K*- . 1 Introduction 

(N 



The possibility to cancel the linear absorption in resonant atomic systems by 

means of electromagnetically induced transparency (EIT) |p| lead in recent 

years to fascinating new developments in nonlinear optics HH] . For example 

Qv' . coherently driven, resonant atomic vapors under conditions of EIT allow for 

complete frequency conversion in distances short enough, such that phase 
matching requirements become irrelevant [0. Furthermore the large nonlin- 
earities of these systems may lead to a new regime of nonlinear quantum op- 

^ . tics on the few-photon level MM with potential applications to single-photon 

C^ I quantum control MM and quantum information processing. 

2^' One particularly interesting nonlinear process based on EIT is the reso- 

nantly enhanced 4-wave mixing in a double- yl system with counter-propagating 
pump modes |£|. It has been shown experimentally ijl^l and theoretically 
[^l]jl^ that this system can show a phase transition to mirrorless oscillations 
for rather low pump powers. Close to the threshold of oscillation an almost 

C^ I perfect suppression of quantum fluctuations of one quadrature amplitude of 

a combination mode of the generated fields occurs [ [l3[[l4[| . Also sufficiently 
above threshold light fields with beat-frequencies tightly locked to the atomic 
Raman-transition and extremely low relative bandwidth are generated |lq| . 
All previous studies of resonantly enhanced 4-wave mixing were done in 
the perturbative regime of small amplitudes of the generated fields. In the 
present paper I want to discuss the case of arbitrary amplitudes. Using a 
simplified open-system model I will derive stationary propagation equations 
for the field amplitudes and present analytic solutions of these equations. 
It will be shown that in an ideal case complete conversion can be achieved 
within a relatively small interaction length. 
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2 Model and Atomic Polarizations 

I here consider the propagation of four electromagnetic waves in a medium 
consisting of double- yl atoms (see Fig.l). These waves include two counter- 
propagating driving fields with equal frequencies fd and Rabi-frequencies Qi 
and ]?2 , and two probe fields (anti-Stokes and Stokes) described by the com- 
plex Rabi-frequencies Ei and E2, with carrier frequencies 1/1=^^ + wq and 
V2 = i^d — 'jJo, where wq — Wbi — Wb2 is the ground-state frequency splitting. 
The fields interact via the long-lived coherence on the dipole-forbidden tran- 
sition between the metastable ground states 61 and &2- We assume that the 
driving field f2i is in resonance with the 62 -^ ai transition, whereas the sec- 
ond driving field J72 has a detuning A^ | J72 1 from the 61 -^ 02 transition. In 
this case linear losses of the fields due to single-photon absorption processes 
are minimized. 




Fig. 1. Atoms in double A configuration inter- 
acting with two driving fields (^1,2) and two 
generated fields (-Bi,2) 



Due to coherent Raman-scattering the pump fields generate counter- 
propagating anti-Stokes and Stokes fields. For a sufficiently large density- 
length product of the medium and for a certain pump field intensity, the 
system shows a phase-transition to self-oscillations ||T^ . The feedback mech- 
anism required for an oscillation is provided here by the gain medium: A 
Stokes photon spontaneously generated on the 02 -^ 62 transition propagates 
in the —z direction and stimulates the generation of an anti-Stokes photon. 
This anti-Stokes photon has a different frequency but a fixed relative phase 
and propagates in the +z direction. It stimulates the generation of another 
Stokes photon upstream. The second Stokes photon will be in phase with the 
first one, provided that the system is approximatly phase matched and that 
there has been no decay of the Raman coherence. The phase-locked emission 
of the second Stokes photon then closes the feedback loop. We have shown in 
[15I that phase-matching enforces a strong pulling of the beat-note of gener- 
ated and pump fields to the atomic Raman transition. I will therefore assume 
here that both A systems are in perfect two-photon resonance. 

In order to calculate the medium response to the fields, one would have to 
solve the atomic density matrix equations to all orders in all fields taking into 
account all relaxation rates. Although this is in principle possible it leads to 



Resonant 4WM: analytic solutions 3 

extremely involved expressions. Instead I here use a simplified open-system 
model which allows to derive rather compact expressions for the atomic sus- 
ceptibilities. 

Since the effects of spontaneous emission are negligible in the present 
system, we may model all relaxations out of the excited states ai and 02 
by rates 7 out of the system. In thermal equilibrium, i.e. in the absence 
of all fields, both lower states 61 and 62 are equally populated. I therefore 
assume - within the open-system approach - that the atoms are pumped 
into states 5i or 62 with 50% probability respectively. The corresponding 
rate is denoted as r and will later be determined by the requirement that 
the total probability to find an atom in any of the states is unity. The finite 
lifetime of the lower-level coherence will here be described by a decay out 
of all states with rate 79. Thus the open-system model corresponds to the 
experimentally relevant situation of an atomic beam or a finite-temperature 
vapor with time-of-flight broadening. In this case the system can be described 
by generalized Schrodinger-equations for field amplitudes instead of density- 
matrix equations. 

The interaction Hamiltonian of an atom at position z with the fields can 
be written in the form 

+Ei{z)e-'''''\ai){h\+E2{z)e~'''-'\a2){b2\+Sidi]. (1) 

If we denote the state vector of the atom as 

\>P) = aie-'""^*|ai) + a2e-'(""^-^)*|a2> + feie-'""*|&i) + &2e-''^-*|62), (2) 

where hv^^ are the energies of the corresponding states, we find the following 
equations of motion of the slowly-varying state amplitudes for an atom at 
position z 

di = -{jQ + j)ai+if2ib2 + iEibi, (3) 

02 = - (70 + 7 + i^) 0.2 + i^2bi + ii52&2, (4) 

bi = ri - ^ah + '^^2^2 + ^E^ai, (5) 

^2 = T'2 - 7o&2 +i^rai +i-B2«2- (6) 

Here I have introduced the rates ri and r2 to distinguish the cases of pumping 
into 61 (ri = r, r2 = 0) and into 62 (^1 = 0, r2 ^ r). Note that simultaneously 
setting ri = r2 = r corresponds to a coherent preparation of the atoms 
in a 50-50 superposition of 61 and &2- In order to describe an incoherent 
preparation in these levels one has to consider the two cases separately and 
add the density matrix elements following from both cases. 

Solving (§-g) in steady state for the case of injection into &i, i.e. for ri = r 
and r2 = one finds 

(1) _ . f2in2E^ - Ei\E2\'' 



\QiQ2-EiE2\ 



(7) 
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(1) _ . 1^1 1 ^2 — f^lEiE2 , , 

^^ =-^^^^^72^3^^^' (10) 

where I have used that Z\ 3> 7 S> 70 and have kept only the leading terms. 
Similarly one finds for injection into 62, i.e. for ri = and r2 = r: 

^(2) _ .^fh\n2l_-EiE2^ 
\QlQ2-ElE2 

(2)_ . Qin2El-\Ei\''E2 



" ion 1:1 1:1 12 ' V-^-^y 



|l2ir22-Sl^2P 



(12) 



'^ =-''- \n,n2-E,E2r ^''^ 

^2 -''^\f2,f22-E,E2\'- ^^^> 

Taking into account only the leading order contribution in the above expres- 
sions is essentially equivalent to assuming an infinitely long lived ground-state 
coherence between 61 and 62- In vapor cells with coated walls or by using 
buffer gases, lifetimes of Hyperfine coherences in alkali vapors in the millisec- 
ond regime are possible. Hence neglecting contributions from finite values 
of 7o seems justified. However, in this case also linear absorption losses are 
neglected. As a consequence the threshold condition becomes independent on 
the pump intensity and an arbitrarily small flux of pump photons is sufficient 
to maintain oscillations p5[ . If on the other hand a small but finite ground- 
state dephasing rate is taken into account, the threshold condition does de- 
pend on the pump intensity leading to a lower limit of the pump-photon flux. 
In the present paper I am interested only in the analytic behavior of the fields 
in the ideal limit and therefore the small but finite linear losses associated 
with the ground-state dephasing will be ignored. 

The pump rate r can be determined from the normalization condition 

Ef^e^'i + e?J = 1- One finds r = (|f2if22 - E^E2f) l\A[\nxf -I- \Ex\^)\. 
With this one obtains for the non-diagonal density matrix elements Qa h — 

\[2^\^Q,[22E*2~ElE2Qin*2 \Qi\' {\Q2? - \E2?) ^ 



ElE*-\E,\^E,E2f2*2 I |giP(|^2p + |i^2p 

i,\' + \Ei\^) E,E2ni \n^m\n,\^ + \Ei 



f2lf22ElE*-\E,\^E,E2f2*2 , l^iPd^zP + |i?2p) ^ ,^.. 

{\Q,\' + \E,\^)E,E2Qt , \Q,?{\Q,? + \E,\^) ^ 
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Qa2b2 



(1^1 



\E,\^)[2,[22El , \Ei\mf2i\ 



\Ei\ 



E. 



(18) 



The first terms in these expressions describe the nonlinear coupling between 
the modes and the second ones ac-Stark shift induced changes in the refractive 
indices. It should be noted that there are no imaginary hnear susceptibilities, 
i.e. there is no linear dissipation despite the fact, that J?i and Ei are in 
single-photon resonance. 



3 Stationary field equations and analytic solutions 



In slowly-varying amplitude and phase approximation, the field amplitudes 
obey the following equation of motion 



dz 



ikiEi 



—E; == i/c2^2 

dz ^ ^ 

d 



dz 



ill — ik^Hi + i 



dz 



■ P^fei 

^2heo 

■ P^fcd 
2heo 

■d 



^2heo 



1- 



N Qaibi, 

N gaib2, 

N Q*a2bi^ 



(19) 
(20) 
(21) 
(22) 



where fci, fc2 and fej are the free-space wavenumbers of the generated and 
pump fields, N is the atomic number density and p are the dipole moments 
of the corresponding transitions, which have been assumed to be equal for 
simplicity. Since the wavenumbers of the fields differ only slightly, one may ap- 
proximate the coupling parameter in all equations by k = p^k^N I2Ti£q. Intro- 
ducing field amplitudes which are slowly varying in space, Ei — Ei e'^^^^, E2 = 
E2 e"'*^^^, ill = ill e''^''^ and il2 = il2 e^''^''^ one eventually arrives at 

I ill I i^l i^2 E2 El E2 i^l i i 2 



-—El = -IK- 

dz 






Ak 



{\ili 



Ai\fli\^ + \Ei\^Y 

^ Ji2i\mi22\'-\E2\') 

" A{\ni\^ + \Ei\^f 

\Ei\^) f2*iil^Ei 



El, 



(23) 



A{\ili\^- 
\Ei\H\ni\'- 



\Ei\^y 

\Ei?) 



Ai\f2i\ 



EX 



dz 



ill 



il-y 1/2 E-y E2 — I El I E\ E2 J i 2 



+!/€- 



A{\f2i\ 
Ai\ili\^- 



\Ei\^) 

\E2n 



\Ei\ 



Hi 



(24) 



(25) 
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+!« ' ' ^' ' ^-^Q^, 26 

where I have dropped the tildes again for notational simplicity, and Ak = 
^2 — fci is the free-space phase mismatch. Expanding these expressions into 
third order of the generated fields Ex and E^ reproduces the equations of 
[fsl. Equations (ES-pq) together with the boundary-conditions 



^i(0)=0, E2{L)^Q, f2i(0) = r2io, and l^aW = ^20, (27) 

where L is the length of the interaction region and Qm and i?20 are the 
given input amplitudes, form a nonlinear boundary- value problem. One easily 
verifies that the set of differential equations has always the trivial solution 
El = E2 = 0, and Qi{z) = ]7io and i72(-z) = ^20- 

As has been discussed in detail in |15[ , the phase mismatch is easily com- 
pensated in an optically dense vapor by a small detuning from the two-photon 
resonance. Oscillation occurs at frequencies such that the phase-matching 
condition is automatically fulfilled. I therefore set this term equal to zero in 
the following. 

Constants of Motion: The field equations have the following constants of 
motion. From the energy-momentum conservation follow the Manley-Rowe 
relations 



_d_ 

dl 
_d_ 

dl 



|r2i|2 + |i?i|2)=0, (28) 

|^2p + |i?2p)=0, (29) 



which state that each photon taken out of the pump fields J7i or f22 is put 
into the anti-Stokcs and Stokes fields Ei and E2 respectively. Furthermore 
one finds that the total intensity of the pump field is constant in space 

^(|f2i|2-t-|/?2|')=0. (30) 

The same is true for the generated fields, which however follow already from 
the above constants of motion. 

'^ ''Eil^ + lE^n^O. (31) 



d. 



z 



Without the phase terms in (g3|-£6|) , which represent contributions due to ac- 
Stark shifts, also the quartic expression Re [/2i/22^i^2] would be a constant 
of motion. In fact the boundary conditions for the generated fields imply that 
Re [f2if22EiE2] = 0. It wiU be shown later on that Re [f2if22ElE2] is in any 
case to a very good approximation a constant of motion. 
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Amplitude-Phase Equations: It is convenient to rewrite the field equa- 
tions in terms of amplitudes and phases. Introducing En = 6„ e~"^" and 



(32) 
(33) 
(34) 
(35) 



i7„ = a„ e "^ 


> (n = l,2) 


one obtains 


d 

d-z'' = 


K O1O262 

Aal + ej 


sini/), 


d 

dz 


K O1O261 

Aaj + el 


sini/), 


d 

-— Ol = 

dz 


K O26162 

Aaf + ej 


sin 7/), 


d 

— 02 = 

dz 


K O16162 

Z\ of + e? 


sini/), 


where tp — cpi 


+ 02-V'l- 


- 1/'2 is the re 


the equation 







d K 

d^^ = Z 



010262(01 - 6^) O1O261 026162(01-6^) 

ei(o2+62)2 

016162 



62(0? +6?) 01(0^ +6^) 



2^2 



02(01 +ef) 



cos ■0 + -r 



2e?o^ 



(of 



o2A2 



(36) 



Solution for Equal Input Intensities: Let me now consider the case of 
equal input intensities of both pump fields, i.e. ai(0) = oio = 020 = cl2{L). 
Making use of the constants of motion one can write 



6i(z) = 6 sin??(z), Oi(z) — Y OiQ — 6^ sin^ '^(•z): 

62(z) — 6 cosi3{z), 02(2) = \/ Oio ~ 6^ cos^ i9{z), 



(37) 
(38) 



with the output amplitude of the generated fields e and the mixing angle i}{z) 
as the only remaining variables. The boundary conditions are now 'i?(0) = 
and ^(i) = 7r/2, if 6 7^ 0, i.e. for the non-trivial solutions. 

Substituting the above expressions into (|32| ) yields the nonlinear equation 



dz A 



1 - £2 + ^ sin^ [2^(z)] 



1/2 



sm^{z), 



(39) 



where e = e/aio. In order to solve ( p9D one can in principle introduce a 
nonlinear stretch of the spatial coordinate according to 



^{z) — / dz' sin ip{z'), and d£^ — sin ip{z) dz 



(40) 



which removes the term sinip^z) on the r.h.s. of (p9|). I will show later on, 
however, that to a very good approximation sin')jj{z) = 1. Thus ^ = z and 
sinip{z) = 1 is used in the following. 
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Integrating ( |39D from z — to z 
ized output amplitude e — e/aio: 

=■4 



K 



4(e2 



1 



f^ 



L leads to an equation for the normal- 



(41) 



where K is the complete elliptic integral of the first kind |T6| . One easily ver- 
ifies that ([4l| ) has only a real- valued solution e, if kL/A > n/2, which is the 
threshold condition for mirrorless oscillations [[lOyi2[] . For smaller values of 
kL/A the equations of motion have only the trivial solution. Figure 2 shows 



\a 



lOl 



0.6 

0.4 
0.2 




12 3 4 5 

KL/A {Till) 

Fig. 2. Output intensity of generated fields E = Ei{L) = E2{0) normalized to 
input intensity of pump fields as function of effective interaction length, i7i(0) = 

n2(L) = nio 



the output intensity of the generated fields normalized to the input inten- 
sity of the pump fields as a function of the effective density length product 
kL. One clearly recognizes that for a sufficiently large product kL complete 
conversion can be achieved. 

The spatial behavior of the field strength inside the vapor cell can be 
obtained from incomplete elliptical integrals following from (p9[). Figure 3 
shows the field amplitudes inside the medium for s = 0.2, i.e. just above 
threshold and for e = 0.98 i.e. for almost complete conversion. 

Not to far above threshold, the square root in (p9|) can be expanded and 
one recovers the third-order solution obtained in flSll: 



^(z) 



KZ 

~A 



1 



1 



with 



V2 



1 



TT A 

2 kZ 



1/2 



for 



kL tt 
~A-2' 



(42) 



(43) 



In order to verify the approximation sin ^'(2:) = 1,1 have numerically inte- 
grated the differential equation (|3^) with the above solutions. Figure 4 shows 
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Fig. 3. Field amplitudes inside interaction region for small conversion e = E/ilio 
0.2 (left) and large conversion e = E/ilio = 0.98 (right) 



the comparison between the nonlinear coordinate ^(z) and z for the case 
e — 0.98. One recognizes that ^ deviates from z by at most 1%. For smaller 
conversions an even smaller difference shows up. Therefore the approxima- 
tion sinV' = 1 is very well justified. This also implies that Re [i7ii72-Ei-E2] i^ 
to a very good approximation a constant of motion. 




Fig. 4. Effective interaction distance ^ versus physical interaction distance z for 
large conversion (e = E/Oiq = 0.98). Dotted line corresponds to ^ = z 



4 Summary 



In the present paper all-order atomic susceptibilities for resonantly enhanced 
4- wave mixing are presented and field equations derived. The coupled nonlin- 
ear differential equations are solved analytically for the case of infinitely long- 
lived ground-state coherences and under the assumption of negligible phase 
changes due to ac-Stark shifts. Below a certain critical value of the density- 
length product only the trivial solution exists, where the generated Stokes 
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and anti-Stokes components have vanishing aniphtude. Above the thresh- 
old to mirrorless oscillations the photon conversion efficiency increases very 
rapidly and at a density-length product of about 3 times the threshold value, 
95% conversion is achieved. 
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